Elimination for generic sparse polynomial systems* 



Maria Isabel Herrero", Gabriela Jeronimo"'^'*, Juan Sabia^'* 

, jj Dcpartamciito de Matematica, Facultad de Ciencias Exactas y Naturales, 

' Universidad de Buenos Aires, Ciudad Universitaria, (1428) Buenos Aires, Argentina 

, t Departamento de Ciencias Exactas, Ciclo Basico Comiin, 

H ■ Universidad de Buenos Aires, Ciudad Universitaria, (1428) Buenos Aires, Argentina 

a • 

' ' o IMAS, CONICET, Argentina 

March 4, 2013 

O 
< 

(— I ■ Abstract 
> ' 

' We present a new probabilistic symbolic algorithm that, given a variety defined in 

an n-dimensional affine space by a generic sparse system with fixed supports, computes 
the Zariski closure of its projection to an ^-dimensional coordinate affine space with 
£ < n. The complexity of the algorithm depends polynomially on combinatorial 
^ invariants associated to the supports. 

CN ! 1 Introduction 

p 

! Let A = {Ai, . . . ,Ar) be a family of finite subsets of (Z>o)" and f = (/i,...,/^) a 

' system of polynomials in Q[Xi, . . . , Xn] supported on A. If V{i) C denotes the afhne 

variety of the common zeros of the polynomials in f and, for a given i < n, it : C" — )• 
is the projection 'ir{xi, . . . , = (xi, . . . , Xi), we consider the problem of computing 



X 



algorithmically a description of the Zariski closure 7r{V{f)) C within a complexity 
depending on combinatorial invariants associated to the input supports. 
. The computation of (Zariski closures of) linear projections of varieties is the basic task 

in elimination theory. A more general formulation of this problem is algorithmic quantifier 
elimination over algebraically closed fields (see, for instance, j40j . |20| . [9], |15] . |32j for 
algorithms with complexities depending on the number and degrees of the polynomials 
and the number of variables involved). Particular instances of the computation of Zariski 
closures of projections are the computation of Chow forms (see, for instance, [7], |18j . 
[25]), classical resultants (see |13] and the references therein) and sparse resultants (see, 
for example, [39], [8], [12], [27]). 

The foundations of the study of sparse polynomial systems can be traced back to [3], 
|28| and [29] . which prove that the number of isolated roots in (C*)" of a system of n 
polynomial equations in n unknowns is bounded by the mixed volume of the family of 
their supports. 
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This result led to the construction of algorithms for the computation of these isolated 
solutions, particularly designed to deal with sparse systems (see, for instance, [12], [23], 
|34j . [26]), which run faster than the known general procedures solving this task. Later on, 
upper bounds for the number of isolated affine solutions of sparse systems were obtained 
and new efficient algorithms to compute them were designed (see, for instance, [33], |30j . 
|35j . |24| . |14j . |16| . |21]). More recently, positive dimensional components of affine varieties 
defined by sparse systems were also considered: in [41j . certificates for the existence of 
1-dimensional components were given and, in [1] and [2], under certain assumptions on 
the equations, algorithmic methods to describe Puiseux series expansions of curves and 
arbitrary positive dimensional components, respectively, were presented. Also, in [22], 
an upper bound in terms of mixed volumes for the degree of the affine variety defined 
by a sparse polynomial system of n equations in n unknowns was proved and algorithms 
for characterizing the equidimensional decomposition of affine varieties defined by sparse 
systems were designed. The good performance of the most efficient algorithms dealing with 
sparse systems relies on the use of polyhedral deformations introduced in [23], because 
these deformations essentially preserve the monomial structure of the system involved. 

In this paper, we present a probabilistic symbolic algorithm which computes the Zariski 
closure 7r(F(f)) C for a sparse system f with fixed supports and generic coefficients. 
We use the decomposition of the variety V{i) into equidimensional subvarieties contained 
in coordinate subspaces proved in [22] to reduce the problem to the case of a variety such 
that each of its components intersects the torus. In this case we compute a geometric 
resolution of the variety with respect to a suitably chosen set of free variables and, from 
this resolution, we show how to obtain a geometric resolution of the Zariski closure of 
the required projection. The complexity of our algorithm is polynomial in combinatorial 
invariants associated to the supports of the input polynomials. Our main result is the 
following: 

Theorem 1 Let f = (/i, . . . , fr) be a system of polynomials in Q[Xi, . . . , Xn] with generic 
coefficients supported on a family A = {Ai, . . . ,Ar) of finite subsets of (Z>o)"' such that 
dimC^j^j Aj) > if J for every J C {!,... ,r}, and let V*{f) C C" be the Zariski closure 
of {x G (C*)'" I fj{x) = for every 1 < j < r}. There is a probabilistic algorithm that 
computes a geometric resolution of 7r{V*{i)), where vr : — )■ is the projection to the 
first i coordinates, within 0{r'D{n'^N \og{d){T>'^ + £)-[- V^)) operations in Q up to loga- 
rithmic factors, where V = My(^, A^""'')), N = I^i<j<r#-4j) d := maxj{deg(/j)} and 
£ = My(((^j)j^/i), A'^""^'"^-^)). Here, MV denotes mixed volume, A is the vertex 

set of the standard simplex in M" and the superscript indicates the number of repetitions. 

The stated complexity is due to a polyhedral deformation based algorithm to solve 
generic sparse zero-dimensional systems (|26j). an algorithmic Newton- Hensel lifting ([36]) 
and codification of multivariate polynomials as straight-line programs ([6]). We do not 
take into account the cost of mixed volume computations (see Section [2. 2p . Furthermore, 
we have ignored terms depending on the size of certain combinatorial objects associated 
to polyhedral deformations involved; for a more precise complexity estimate, see Theorem 

m 

The paper is organized as follows. In Section [21 the notation and the basic theoretical 
and algorithmic notions used throughout are introduced. Section [3] is devoted to proving 
the main theoretical results on which our algorithms rely. Finally, Section H] contains the 
descriptions, proof of correctness and complexity estimates of our algorithms and examples 
illustrating how they work. 
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2 Preliminaries 



2.1 Basic definitions and notation 

Let A; be a field of cliaracteristic zero and k be an algebraic closure of k. Given polynomials 
/ii ■ ■ ■ ) /r & k[Xi, . . . ,Xn], we write V{f) = V{fi, . . . , fr) for the affine variety which is 
the set of the common zeros of f = (/i, . . . , fr) in k , and 

y*(f) = v*{h, ...,fr):= y(/i,...,/.)n(r)- 

for the union of the irreducible components of V{fi,... , fr) containing points with all 
their coordinates in A; := k — {0} 

For a variety V C k definable over k, we denote k[V] = k[Xi, . . . , Xn]/I(V) its 
coordinate ring (where I{V) C k[Xi, . . . , Xn] is the ideal of the polynomials vanishing 
identically on V). V is irreducible, we write k{V) for the fraction field of k[V]. 

To describe zero-dimensional affine varieties we use the notion of a geometric resolution 
(see, for instance, |19j and the references therein): Let V = ■ ■ ■ C be a 

zero-dimensional variety defined by polynomials in k[Xi, . . . ,Xn]. Given a linear form 
A = AiXi + • • • + XnXn in klXi, . . . such that X{^^'^) / X{^^i^) if i / j, the following 
polynomials completely characterize V: 

• the minimal polynomial qx = ni<j<D 

(Y - A(e(*))) G k[Y] of A over the variety V 

(where y is a new variable), 

• polynomials vi,...,Vn € k[Y] with deg(t;j) < D for every 1 < j < n satisfying 
V = {{viiri), Vniv)) €C^\v^k, qx{v) = 0}. 

The family of univariate polynomials {qx,vi, . . . ,Vn) G k[Y]^~^^ is called the geometric 
resolution of V (or a geometric resolution of k[V]) associated with the linear form A. 

The notion of geometric resolution can be extended to any equidimensional vari- 
ety: Let V C k be an equidimensional variety of dimension t defined by polynomi- 
als in k[Xi, . . . ,Xn]- Assume that, for each irreducible component W of V, the iden- 
tity I{W) n k[Xi,...,Xt] = {0} holds. By considering k{Xi,...,Xt) (g) k[V], we are 
in a zero-dimensional situation, and we call a geometric resolution of V with free vari- 
ables Xi, . . . ,Xt to a geometric resolution {qx,vt+i, . . . , Vn) G k{Xi, . . . , Xj)[y]'"~*+-^ of 
k{Xi, . . . , Xt)(8>/c[y] associated to a linear form A G k[Xt+i, . . . , Xn]- If G k[Xi, . . . , X^, y] 
is obtained from qx by clearing denominators, a geometric resolution of V gives a bira- 
tional map between the hypersurface {{xi, . . . ,xt,y) G AT^ \ qx{xi, . . . , xt,y) = 0} and 
V. 

When dealing with varieties defined by sparse polynomial systems, an important com- 
binatorial invariant associated to the system is the mixed volume of their supports. For 
a family A= {Ai,..., An) of n finite subsets of (Z>o)"', MVn{A) (or MV{A) if n is clear 
from the context) denotes the n-dimensional mixed volume of the convex hulls in of 
Ai, . . . ,An (see, for instance, [lOl Chapter 7] for the definition and basic properties). In 
this context, we write A for the vertex set {0, ei, . . . , e„} of the standard simplex of M", 
which is the support of an affine linear form with nonzero coefficients, and A^*-* for the 
family of t copies of A. 

2.2 Algorithms and codification 

Although we work with polynomials, our algorithms only deal with elements in a base 
field k. The notion of complexity of an algorithm we consider is the number of operations 
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and comparisons in k it performs. We will encode multivariate polynomials in different 
ways: 



• in sparse form, that is, by means of the list of pairs (a, Ca) where a runs over a 
fixed set including the exponents of the monomials appearing in the polynomial 
with nonzero coefficients and Ca is the corresponding coefficient, 

• in the standard dense form, which encodes a polynomial of degree bounded by d 
as the vector of the coefficients of all the monomials of degree at most d including 
zeroes (we use this encoding only for univariate polynomials), 

• in the straight-line program (sip for short) encoding. A straight-line program is an 
algorithm without branchings which allows the evaluation of the polynomial at a 
generic value (for a precise definition and properties of sip's, see [6]). 

In our complexity estimates, we use the usual O notation: for f,g: Z>o — > M, f(d) = 
0{g{d)) if \f{d)\ < c\g{d)\ for a positive constant c. We also use the notation M{d) = 
(ilog^((i) log(log((i)), where log denotes logarithm to base 2. We recall that multipoint 
evaluation and interpolation of univariate polynomials of degree d with coefficients in a 
commutative ring R of characteristic can be performed with 0{M(d)) operations and 
that multiplication and division with remainder of such polynomials can be done with 
0{M{d)/ log{d)) arithmetic operations in R. We write Q for the exponent (fi < 4) in the 
complexity 0{d^) of operations (determinant and adjoint computations) on dxd matrices 
with entries in a commutative ring i? ([!]). 

Our algorithms are probabilistic in the sense that they make random choices of points 
(that we consider cost-free in our complexity estimates) which lead to a correct compu- 
tation provided the points lie outside certain proper Zariski closed sets of suitable affine 
spaces. Using the Scwhartz-Zippel lemma (|371l43j). the error probability of our algorithms 
can be controlled by making these random choices within sufficiently large sets of inte- 
ger numbers whose size depend on the degrees of the polynomials defining the previously 
mentioned Zariski closed sets. 

Different steps of our algorithms rely on mixed volume computations. To perform 
them we may apply the dynamic enumeration procedure from [31], which proved to be 
efficient even for large systems. However, there are no explicit complexity bounds for this 
procedure; for this reason we do not include its cost in our complexity estimates. 

3 Theoretical results 

Let n,r be positive integers and A = {Ai, . . . ,Ar) a family of r finite subsets of (Z>o)'^. 
Consider r polynomials in n variables X = {Xi, . . . , Xn) supported on A with indetermi- 
nate coefficients: for every 1 < j < r, let 



where Cj := {Cj^a)a<^Aj is a set of Nj := #Aj new indeterminates. Let f := (/i, . . . , fr) 
and K := Q(C7i, '. . . , C^). 




(1) 



aeAj 
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3.1 Reduction to the toric case 

Given the family f = (/i, . . . , /r) of polynomials in K[Xi^ . . . , Xn] supported on A intro- 
duced before, the equidimensional decomposition of the affine variety V{{) C K depends 
only on the combinatorial structure of A (see j22j). Moreover, its equidimensional com- 
ponents can be defined from certain smaller polynomial families f/ in fewer variables 
associated to f . This decomposition, which we explain below, enables us to reduce the 
problem of computing a projection of the affine variety V{i) to the computation of finitely 
many projections of varieties having a non-empty intersection with the torus. 

For I C {!••., n}, given a polynomial / € K[Xi, . . . , Xn], we write // to denote the 
polynomial in K[{Xi)i^j] obtained from / by specializing Xi = for every i ^ I. We 
define J/ = {j € {1, . . . , r} | 3a S : Oj = Vi G /} and f/ = {{fj)i)jeJi (note that Jj is 
the set of indices of those polynomials fj that do not vanish identically when specializing 
Xi = for every i £ I). Let 

r = {/ C {1, . . . , n} I V J C J/, dim(^ A'j) > # J; V7 C /, # + #7 > # J/ + #/}, 

where, for every I < j < r, Aj C (Z>o)"~*^ denotes the support of {fj)i- Finally, let 

: K""*^ ^ be the map that inserts zeros in the coordinates indexed by I. Using 
the previous notation, according to [22\ Theorem 7], we have that 

v{{) = u Mv*m- 

/er 

Let vr : — > be the projection 7r(xi, . . . , Xn) = {xi, . . . , x^). For / C {1, . . . , n}, 
let ii = £ — #(/ n {1, . . . TTij : K"' — > K^' the projection to the first coordinates 
and </?/n{i,.../} • ~^ the map that inserts zeros in the coordinates indexed by 
In {1, . . . ,£}. As TT{ipi{V*{{i))) = V?7n{i,...,^}(^^/(^*(f/))) holds for every /, we conclude 
that 

lev 

Therefore, the computation of 7r(y(f)) amounts to obtaining the Zariski closures of 

the projections by vr^^ of the afhne varieties V*{ii) = {x G (i('*)"~#^ | f/(x) = 0}. For this 
reason, in the sequel we will deal with this problem. 

3.2 The vanishing ideal of the projection 

We consider a generic polynomial system f := (/i, . . . , fr) supported on A := {Ai, . . . , Ar) 
with indeterminate coefficients Cj := {Cj^a)aeAj for every 1 < j < r as in equation ([T]). 
We assume that r < n and that, for every J C {!,..., r}, dim{J2j^j Aj) > #J (or, 

equivalently, My(^, A^'""'')) > 0) so that V*{i) is not empty. 

Our aim is to obtain a geometric resolution of the Zariski closure TT{V*{i)), where 
vr : {£ < n) is the projection defined by 7r(rEi, . . . , x„) = (xi, . . . , x^). We begin 

by proving some basic results on the vanishing ideals of the varieties V*{i) and 7r(y*(f)). 

Lemma 2 Under the previous assumptions, (/i, . . . , fr) ■ {X\ . . . is a prime ideal 

in K[Xi, . . . , Xn] of dimension n - r. Moreover, I{V*{i)) = (/i, . . . , fr) ■ {Xi . . . X„)°°. 
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Proof: For every 1 < j < r, fix ajo G Aj and let A'j := Aj — {ajo} and Cj := Cj — {Cj^ajo}- 
We denote C := (Ci, . . . , Cr) and C := {C[, . . . , C^). Consider the ring morphism ip : 
Q[C,X] Q[C',X]xi...x„ (where Q[C", is the localization of Q[C',X] at the 

multiplicative set | m £ Z>o}) defined by ipiXi) = for 1 < i < n, 

V'lC,,,) = for 1 < j < r, a e A'^, and ^(C,-„^J = -X-^^o^ ^ C7,-,X'^) for 1 < j < r. 

We claim that ker(V') = (/i, . . . , /r) : (^i • • • Xn)°°- It is clear that for a polynomial 
5 G (/i, . . . , /r) : (^1 . . . ^n)°°, ''/'(S') = 0. Assume now that for a polynomial g € Q[C, X] 
we have that ipig) = 0. Let Co := (Ci^aio) • • • ' C'r^aro) ^^"^ := (C'lOi • • • ) C'ro) be new 

variables. By Taylor expansion, g{Co,C',X) = g{Co,C',X) + I]i<j<,,(CjO - Cj,ajJ • 
for certain polynomials Gj G Q[do,C,X]. Specializing djo = Yl Cj aX"-) for 

every 1 < j < r, it follows that ^(5) = g{Co,C',X) - Ei<i<r ^""'"/j ' with Gj G 
Q[C,X]xi...x„; therefore, g{C,X) = Ei<i<r ' G,- Multiplying by (Xi . . 

for a sufficiently large m, we conclude that {Xi . . . Xn)"^g{C, X) G {fi, . . . , fr). 

Then, (/i, • • • , /r) : (^1 • • • Xn)°° is a prime ideal of Q[C, X]. The first statement of 
the lemma follows by localizing at Q[C] — {0}. 

In order to prove the second part of the lemma, consider first a polynomial g G 
(/i, ...,/,): (Xi .. . Xnr, and let m G Z>o be such that (Xi . . . X„)'"5_G (/i, . . . , /,). 
Then, (Xi . . .Xn)'^g vanishes over V{i) and so, (7 vanishes over y(f) n (iv: therefore, 
g G I(y*(f)). Conversely, if 5 G /(y*(f)), then it vanishes over each irreducible compo- 
nent of V{i) intersecting {K )" properly. Then, {Xi . . . Xn)g vanishes over V{fi, . . . , fr), 
which implies that there exists m G Z>o with [Xi . . . Xn)"^g"^ G (/i, . . . , /r). There- 
fore, g'^ G (/i,...,/r) : (Xi . . . Xn)°° and, since this is a prime ideal, it follows that 
5G (/i,...,/.) : (Xi...X„)-. □ 

Corollary 3 The affine variety V*{i) C is an irreducible K -variety of dimension 
n — r. 



Taking into account that for any K-vaviety V C K , the identity I{it{V)) = I{V) fl 
K[Xi, . . . ,Xi] holds. Lemma [2] also enables us to characterize the vanishing ideal of the 
projection we want to compute: 

Corollary 4 With the previous assumptions and notation, 

/(vr(F*(f))) = ((/i, ...,/,): (Xi .. . XnD n K[X,, . . .,X,]. 



Let t := dim(7r(l^*(f))). Without loss of generality, by renaming variables, we may 
assume that {Xi, . . . ,Xt} C {Xi, . . . ,Xi} is a transcendence basis of K{-ir{V*{i))) over 
K and {Xt+r+i, • • • , X„,} are such that {Xi, . . . , Xt, Xt^r+i, ■ ■ ■ , Xn} is a transcendence 
basis of K{V*{{)) over K. The following proposition generically allows us to deal with 
projections with 0-dimensional generic fibers. 

Proposition 5 There is a Zariski dense open set O C 'K^ * ^ such that, for every b G 
K^-'-'-nO, the identity /(7r(y*(f))) = {{f,{Xi, . . . , Xt+r,b), . . . , fr{Xi, . . . , Xt+r,b)) : 
(Xi . . . Xt+r)°°) n K[Xi, ...,Xi] holds. 
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Proof: Let X := {Xi, . . . , Xt+r)- From Corollary HI it is clear that /(7r(y*(f)) C 
{{MX, h),..., fr{X, b)) : (Xi . . . Xt+rD n K[Xi, ...,Xe] for every b = {bt+r+i, ...,bn) 
such that bi ^ for every t + r+ l<i<n. 

For the converse inclusion, first note that, for every 1 < j < r, 

K[Xi, . . . ...,/,■): (Xi .. . X^r ^ K[Xi, . . .,Xr,\x,...x„)/{fu • • • , fj) 

~ K[Y,Xi, . . .,Xn]/{YXi ...Xn-lJu..., fj). 

As in the first part of the proof of Lemma [21 we have that (/i, . . . , fj) : {Xi . . . Xn)°° is a 
prime ideal of dimension n — j for every 1 < j < r; therefore, YXi . . . X„ — 1, /i, . . . , is a 
reduced regular sequence in ^[1", Xi, . . . , Moreover, the set {Xi, . . . , Xt, Xf+r+i, ■ ■ ■ , ^n} 
is algebraically independent modulo {YXi . . . Xn — 1, /i, . . . , /r). Then, by [TTl Corollary 
17 and Theorem 19], there exists a iC-definable Zariski dense open set O C ^ con- 

taining {xi / V 1 < z < i + r} such that for every b £ O D i^""*"'', Cfe YXi . . . Xt+r - 
1, fi{X, b), . . . , fr{X, b) (where Cb := bt+r+i ■ ■ ■ bn) is a reduced regular sequence in K[X] 
and {Xi, . . . , Xt} is algebraically independent modulo each of the associated primes of the 
ideal this regular sequence generates. By noticing that 

6), . . . , fr{X, b)) :{X,... Xt+rr ^ 

~ K[X]/{ci, YXi... Xt+r - 1, fi{X, 6), . . . , fr{X, b)), 

we conclude that {fi{X, b), . . . , fr{X, b)) : {Xi . . . Xt+r)°° is a radical equidimensional ideal 
of dimension t and {Xi, . . . ,Xt} is algebraically independent modulo each of its associated 
primes. Then, the same holds for the ideal {{fi{X,b), . . . ,fr{X,b) ) : {Xi . . . Xt+r)°°) n 
K[Xi, . . . , Xi\. As this ideal includes the t-dimensional prime ideal I {'ir{V *{¥))), the equal- 
ity in the statement of the proposition holds. □ 

3.3 Free variables 

The following result gives a combinatorial condition for the algebraic independence of a 
subset of variables modulo the vanishing ideal of V*{f) that we will use in our algorithm 
to compute a suitable transcendence basis of the rational fraction field of this variety. 

Lemma 6 Let /i, . . . , € K[Xi, . . . be sparse polynomials supported on a family 

A = {Ai, . . . , Ar) of finite subsets of (Z>o)" as introduced in (Cp, and let I = (/i, . . . , fr) : 
{Xi . . . Xn)°° C K[Xi, . . . , Xn]. Let {ei, . . . , 6^} be the canonical basis of Q". 

Then, the set {X,^ , . . . , Xi^, } C {Xi , . . . , Xn} for 1 < ii < ■ ■ ■ < i^ < n is algebraically 
independent modulo I if and only if MV{Ai, ■ ■ ■ , Ar, {0, e,,^ },..., {0, Cj^. }, A^'""'""'^)) > 0. 

Proof: First, assume that MV{Ai, • • • , A, {0, J, . . . , {0, e^^ }, A^""''"*^)) > 0. Let Q be 
a non-zero polynomial in the coefficients of a system supported on {A, {0, Cj^ }, . . . , {0, ej^. }, 
^(n-r-fc)^ such that for each coefficient vector (c, r?) = (ci, . . . ,Cr,r]r+i, ■ ■ ■ ,r]n) with co- 
ordinates in Q* and Q{c,r]) ^ 0, the corresponding sparse system has as many isolated 
solutions in (C*)" as the mixed volume. 

Let p € X n K\Xi^, . . . ,Xi^. Without loss of generality, we may assume that p € 
Q[C, Xi^ , A^ife] • We have that there exist m G Z>o and a non-zero polynomial pq G Q[C] 
such that 

Po(C)(Ai . . . A„)XC, A,, , . . . , AiJ = + • • • + 5,/, (2) 
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with gi,...,gr € Qp, Xi, . . . , 

For each (c, r]) with coordinates in Q* such that pq{c)Q{c, rj) ^ 0, considering a solution 
C G (C*)" of the system 

fl{ci,X) =0,..., fr{Cr,X) = 0, rir+l^l + 77r.+i,2Xij = 0, . . . , ?7,.+fc,l + ^r+fc,2-'^ifc = 0, 

^r+fc+1,0 + ^ Tlr+k+l,iXi = 0, . . . , 7?n,0 + ^ r]n,iXi = 
l<i<n l<i<n 

and specializing identity ^ in (c, 77, ^) we obtain that p{c, — 77r+i,i/??r+i,2) • • • , —ilr+k,i/Vr+k,2) = 
0. We conclude that p = 0. 

Assume now that {Xi-^, . . . ,Xi^ } is algebraically independent modulo I = I{V*{i)). 
Let li, . . . , In-r-k be linear forms in the variables Xi, . . . , Xn with coefficients in K* such 
that {Xi^, . . . , Xii_,li, . . . ,ln-r-k} is a transcendence basis of K{V*{{)). Since V*{{) fl 
{K )" is a dense open subset of F*(f), for a generic (Ci, • • • ,Cn-r) S (-fC )'^~'', we have 
that 

y*(f) n {Xi^ = V 1 < J < A:} n {/^(x) = C^+i V l < j < n - r - A:} 

is not empty and consists of finitely many points in (i^*)". These points are the common 
solutions in (K )" of the system 

h{x) = 0,. . .,fr{x) = o,x,, - Ci = 0, . . . - a = 0, 

hiX) — Ck+l = 0, . . . , ln-r~k{X) — Cn-r = 0, 

which is supported on (^i, . . . , Ar, {0, e^^ }, . . . , {0, e^^. }, A^""^"*^)). By Bernstein's Theo- 
rem, we conclude that My(^i, . . . , A, {0, J, . . . , {0, J, A^"-'^"'^)) > 0. □ 

Remark 7 For 1 < ii < • • • < ifc < we have that 

MVniAi, . . .,Ar, {0,e,J, . . . , {0,e,J, A("— = MK_fc(t77(A), . . . ,w{Ar),w{A)^^-'~^^) 

(see f38l Lemma 6]), where the projection to the coordinates indexed 

by {1, . . . ,n} — {ii, . . . ,ik}- This implies that, in order to determine whether a set of k 
variables is algebraically independent or not, it suffices to compute an {n — k)- dimensional 
mixed volume. 

4 Algorithms 

In this section we will present an algorithm to compute the Zariski closure of the projection 
of y*(f) to the first i coordinates, where f is a generic polynomial system with given 
supports. First, we describe some subroutines. 

4.1 Subroutines 

The first subroutine we will use, which follows from Lemma [6l computes a transcendence 
basis of K{V*{{)) containing a transcendence basis of K{7:(y*{i))). 



Algorithm TransBasis 

INPUT: A family A = {Ai, . . . , Ar) of finite subsets of (Z>o)" such that dim{Y,jizj Aj) > 
# J for every J C {1, . . . , r}. 
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1. TB := 



2. k:=l 

3. while i^TB < n — r do 

(a) Compute D := MV{Ai, • • • , A, ({0, e^J) .^^^^ , {0, e^}, A^"-'^"*^^"!)). 

(b) If L» > 0, := TB U {/c}. 

(c) k:=k + l 

OUTPUT: The set = {ii, . . . , in-r} with ii < • • • < such that {Xi^,. . . , Xi^_^} 
is a transcendence basis of K{V*{f)) over and, for each 1 < j < n — r, {Xi^^, . . . ,Xi.} 
is a maximal algebraically independent subset of {Xi, . . . ,Xi^}. 



Without loss of generality, by renaming variables, we may assume that the tran- 
scendence basis of K(y*{{)) obtained by applying algorithm TransBasis is {Xi, . . . ,Xt, 
Xt+r+i, • • • ) Xn} with t < i and i < t + r. Then, {Xi, . . . , Xt} is a transcendence basis 
of K{7r{V*{{))). Let be the polynomial system obtained by evaluating Xt+r+i, ■ ■ ■ ,Xn 
in a generic point b. By Proposition [5l we may obtain 7r(y*(f)) using the system f^. In 
order to do this, we will compute first a geometric resolution of V*{fb) by means of the 
subroutine we introduce below. 

Let A; be a field of characteristic 0. Algorithm ParametricToricGeomRes computes 
a geometric resolution of the variety V*{g) defined from a generic sparse system g := 
{gi, . . . , Qr) in k[Xi, . . . , Xt+r] with given supports S := . . . , Sr), provided that {Xi, . . . , Xt} 
is a set of independent variables for all its components. This subroutine is obtained 
by following the parametric geometric resolution algorithm from |361 Theorem 2] taking 
Xi, . . . ,Xt as the parameters. 



Algorithm ParametricToricGeomRes 

INPUT: A generic sparse system g := {gi, . . . ,gr) in k[Xi, . . . ,Xt+r] with given supports 
S := {Si, . . . ,Sr) such that {Xi, . . . ,Xt} is algebraically independent modulo each asso- 
ciated prime of {gi,... ,gr) : (^i . . . Xt+r)°°. 

1. Choose ^ = (^1, . . . ,S^t) at random with € Z — {0} for every 1 < i < t. 

2. Compute a geometric resolution of the common solutions of g(^, Xf+i, . . . , Xt+r) in 

(ky. 

3. Obtain an sip encoding the polynomials in g. 

4. Apply a symbolic Newton-Hensel lifting (in the parameters Xi, . . . ,Xt) to the geo- 
metric resolution obtained in step 2 with precision 

5. Applying Pade approximation to the output of the previous step, recover numerators 
and denominators in k[Xi, . . . , Xt] for the coefficients of the polynomials in the 
geometric resolution of U*(g). 
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OUTPUT: A geometric resolution of V*{g) with free variables Xi, . . . ,Xt. 



Now we explain how the different steps of the algorithm are carried out and estimate 
the number of operations in k it performs. We will use the notation = X^x<j<r #Sj 
and d = maxi<j<r deg(gj). 

In Step 2, the algorithm computes the sparse encoding of g{^,Xt+i, . . . ,Xt+r) within 
0{r{t + r)N\og{d)) operations. Then, this system is solved using the procedure from [26] 

within 0{r^Nlog{d)MiM)MiMVim{S))){MiMV{w{S)))+M{u^^^Y:i^^^ 
where zu is the projection to the last r coordinates, Wmax is the maximum of the val- 
ues taken by a generic lifting function uj for w{S) := {vj{Si), . . . ,w{Sr)), and A4 := 
max{||/x||}, the maximum ranging over all primitive normal vectors to the mixed cells in 
the fine mixed subdivision of zu{S) induced by uj. 

In Step 3, an sip of length 0{{t + r)Nlog{d)) encoding the polynomials g is obtained 
from their sparse representation. 

The next step is performed modifying the procedure in [36( Section 4.2] in order 
to use straight-line programs for computations with truncated multivariate power se- 
ries. We represent each of these series as the vector of its homogeneous components 
and these components by means of straight-line programs. The required precision is 
2deg(y*(g)) = 21/1^(5, AW) (see [22l Lemma 1]). Therefore, this step is done within 
complexity 0{{r{t + r)Nlog{d) + r^)M{MV{w{S)))MV{S,A'^*'if) and produces an sip of 
the same order encoding the homogeneous components of the coefficients of the output. 

Finally, the Pade approximation to obtain the coefficients of the geometric resolution is 
done by reducing it to a univariate problem and solving it by means of subresultant com- 
putations following [171 Corollaries 5.21 and 6.48]. This step adds 0(rMy(5, AW)^+2) 
operations to the previous complexity and 0(rMy(5, A^*))^"*"^) to the sip length. 

The total number of operations in k performed by the algorithm is of order 

0((r^ + rt)Nlog(d)MiM)M(MV(w(S)))(MV(S,A^^h^ + 

(3) 

+ M(^^ax Ei</.<r MV{w{{Sj)j^h),A))) + rMV{S, A(*))^+2) 

and produces an sip of length 

0{rMV{S,A^^'^)^{{t + r)Nlog{d)+r^)M{MV{w{S))) + MV{S,A^^^f-^)) (4) 

for the coefficients of the output. 

The last step of our main algorithm consists in describing the projection to a coordinate 
subspace of an equidimensional variety of dimension t given by a geometric resolution in 
the case that the projection has the same dimension t. To do this, we apply the subroutine 
described below. 

— t-\-r 

Let V C k be an equidimensional variety of dimension t definable over k and such 
that for each irreducible component W of V, I(W) n k[Xi, . . . , Xt] = {0} holds. Consider 

the projection it : k k where £ > t, tt{xi, . . . ,xt+r) = {xi, ■ ■ ■ ,X£). Note that 

{Xi, . . . ,Xt} are free variables with respect to each irreducible component of 7r{V). Let 
K := k{Xi, . . . ,Xt). Suppose A € k[Xt+i, . . . ,Xt+r] is a primitive element for K (g) k[V] 
and let {qx,wt-\-i, . . . ,wt+r) G ]K[y]^"'"^ be the associated geometric resolution. Let D be 
the dimension of IK (8) k[V] as K- vector space. 
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Let /.f = f-H+iXf+i + • • • + niXi be a primitive element for K (g) A:[7r(y)]. As I{-k{V)) = 
I{V) n k[Xi, . . . , Xi], to find the minimal polynomial of /i with respect to 7r(y), it suffices 
to find a polynomial € IfC[y] of minimal degree such that G IK ® /(V^). Then, 

6 := degy(g^) is the dimension of IK (8) /c[7r(y)] as a IK-vector space and so, for each 
t + 1 < j < £, in order to obtain a polynomial vj such that Xj = Vj{n) in IK (g) A;[7r(y)] it 
suffices to find a linear combination Xj = Yli=o "^jil^^ of {!) A*) • • • > A**^"^} in IK C?) k\y]. To 
do this, we use the basis B\ := {1, A, ... , A^~^} of IK (g) k\y]. 

In order to compute the geometric resolution of vr(y) associated with ^u, we first look for 
the minimal power which is a IK-linear combination of {1, ^, . . . , in IK(8)A;[V]. Since 

Xj = Wj{X) for every t + 1 < j < t + r, we have that /x = fJ-jWj{X) = p^{X), where 

P^^(Y) := Ej=t+i t^jWjiY) and, for every ieN, fi' = p^(A)* = {Pf^iYf mod gA(>'))|y=A- 
Therefore, 5 equals the rank of the D x D matrix whose columns are the coefficients of 
the polynomials {p^{Yy mod q\{Y)) for z = 0, . . . , D — 1. 

Then, we obtain the minimal polynomial qfj,{Y) := + X]i=o 9^,1^* of /.i and the 
polynomials Vj{Y) = YliZo VjiY^, t + 1 < j < i, which form the geometric resolution of 
7r{V), by solving the linear systems obtained by equating the coefficients of the different 
powers of A in the identities 

<5-l 5-1 
Pf^sW = ^{-qfj.,i)Pti'W and Wj{\) = ^Wjip^i(A), t + l<j <£. 

i=0 i=0 

Summarizing, with the previous notation and hypothesis, we have: 



Algorithm GeomResProj 

INPUT: A geometric resolution {qx, wt+i, . . . , wt+r) of V with free variables Xi, . . . , Xt, 
and a linear form fj, = Yl^j^i l^t+jXt+j € k[Xt+\, . . . , X(\ which is a primitive element for 

^(^k[K{V)]. 

1. Setp^o(y) := 1 andp^(y) := Ylh=oiY!'j=t+iH'^j,h)^^ ^ where {wj^, . . .,Wj^D-i) =■ 
Wj is the vector of coefficients of Wj{Y) for j = t + 1, . . . ,t + r. 

2. For i = 2,...,D, compute p^.{Y) := (p^(y) -p^.-iiY) mod qx{Y)) 

3. Compute 6 := rank(p^jO, p^, . . . , p^o-i), where p^i G IK-^^^ denotes the vector of 
coefficients of p^i . 

4. SetP:= (p^o,p^,...,p^.-i) GK^x^. 

5. Solve the linear systems P • q = p^* and P ■ Vj = Wj, t + 1 < j < £, to obtain 
q := {qo, qs-i) and := {vj^, Vj^s-i) 

6. Set q^iY) := Y^ - Ylto Qi^' and Vj{Y) := Ylto VjiY' for every t + 1 < j < 1 

OUTPUT: The geometric resolution {q^,vt+i, . . . ^vi) of the projection it{V) C k^ with 
free variables Xi, . . . ,Xt associated to the linear form /i. 



The correctness of the procedure follows from our previous arguments. Now we esti- 
mate its complexity. Assume that the input polynomials are encoded in dense form as 
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degree D univariate polynomials in k(Xi, . . . , and their coefficients are encoded by 

an sip over k of length L. 

In Step 1, the algorithm computes an sip encoding the coefficients of of length 
boundedby L + 2L'(^-t). 

In order to fulfill Step 2, we first compute recursively the powers y'^ for /i = L>, . . . , 2D- 
2 modulo q\{Y) from the coefficients of qx and then, we obtain an sip of length 0{L+D{1— 
t)+D'^) for the coefficients oip^i for i = 2, . . . , Z) by expanding the product Pi_i{Y)-p^i-i (Y) 
and substituting the powers of Y by their previously computed expressions. 

Step 3 is done probabilistically by choosing a point x = {xi, . . . ,xt) at random, evalu- 
ating the involved rational functions at this point within 0{L + D[£ — t) + D^) operations 
in k and finally computing the rank 5 of (p^o(2;), p^(2;), . . . ,p^D-i(rc)) within 0{D^) op- 
erations in k, with uj < 3 (see [3 Chapter 2, Sec. 2]). 

To solve the linear systems involved in Step 5, the algorithm computes the invertible 
matrix P*P within 0{D5^) operations in k, its adjoint matrix and determinant within 
0{5^) additional operations, and the products adj(P*P)p^« and adj(P*P)wj for t + 1 < 
j <t with 0{6'^{l-t)). 

Adding the previous estimates, we conclude that the algorithm produces an sip of 
length 0{L + D^ + D6{i -t) + 6^) over k within a complexity of the same order. Taking 
into account that 6 < D, it follows that both the number of operations in k performed by 
the algorithm and the length of the output sip are of order 

0{L + + D^{i-t)). (5) 
4.2 An algorithm to find tiie projection 

Here we present a probabilistic algorithm that, from a fixed family A = {Ai, . . . ,Ar) of 
finite sets (Z>o)"' and a fixed family of variables Xi, . . . , Xi, obtains a geometric resolution 

of 7r(y*(f)) C K , where f = (/i, . . . , fr) is defined in equation ([T]) and tt : K K is 
the projection 7r(xi, . . . , Xn) = {xi, . . . , Xi). 



Algorithm iv -Projection 

INPUT: A family ^ = (^i, . . . , A) of finite subsets of (Z>o)'' such that dim{Y,j(.j Aj) > 
#J for every J C {1, . . . , r} and a set of variables {Xi, . . . , Xi} C {^i, . . . , Xn}- 

1. Apply Algorithm TransBasis to the family A. Without loss of generality, we suppose 
the transcendence basis obtained is {Xi, . . . , Xt, Xt+r+i, ■ ■ ■ , Xn} with t < i and 
t + r + l>i. 

2. Choose randomly in Z the entries of a vector b = {bt+r+i, • • • , bn) and of a vector 
(At+i, . . . , Xt+r)- 

3. Obtain the sparse representation of the system of polynomials 

ffo = ■ ■ ■ , Xt+r, b), . . . , fr{Xl, . . . ,Xt+r, b)) in K[Xi, . . . , Xt+r]- 

4. Apply Algorithm ParametricToricGeomRes to the system and the variables Xi, . . . ,Xt 
to obtain the geometric resolution (qx, wt+i, - - - , wt+r) of the variety V*{ib) with free 
variables Xi,...,Xt associated to the linear form A = Xt+iXt+i -|- • • • -|- \t+rXt+r- 

5. Choose randomly in Z the entries of a vector (^t+i, . . . , /i^). 
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6. Apply Algorithm GeomResProj to {qx,Wt+i, . . . , Wt+r) and fi = fit+iXt+i + - • ■+fiiXi. 

OUTPUT: A geometric resolution (g^, vt+i, . . . ,ve) of 7r{V*{i)) C , where vr : 
is the projection to the first coordinates. 



Theorem 8 Given a family A = (^i,...,^^) of finite subsets of (Z>o)" such that 
d\m.{^-^j Aj) > ij^J for every J C {l,...,r} and the projection it : — )■ to the 
first coordinates, Algorithm i^T-Projection is a probabilistic procedure that computes a 
geometric resolution o/7r(y*(f)) for the sparse system f supported on A with indetermi- 
nate coefficients within 

0{{n^ + r^)N\og{d)M{V)E{V^ + M{£)) + rV^+^) 

operations in K , where N = Ylii<j<T T^^i> ^ •= ™a'Xi<j<f,{degx(/j)}, ^ = MV{A, A^""^)), 
8 = X]i<h<r A^""*""^^)), and 3 is a constant measuring the size of certain 

combinatorial objects involved at intermediate computations and associated to the family 
of supports A. 

Proof: Since in our complexity estimates we only take into account the number of op- 
erations in K (and not mixed volume or mixed subdivision computations), to obtain the 
overall complexity of the algorithm it suffices to add the complexities of steps 3, 4 and 6. 
Step 3 can be done within 0{n'^Nlog{d)) operations in K. 

The complexity of Step 4 is the already stated for Algorithm ParametricToricGeomRes 
in ([3]) . Note that for generic 6, the system fb is a generic polynomial system supported on 
S := {Si, . . . , Sr), where Sj C {Z>qY~^^ is the projection of Aj to the first t + r coordinates 
for every 1 < j <r. Moreover, by Bernstein's theorem, My(5, A^*)) < MV{A, A^^~^^), 
MV{w{S)) < My(^, A("-^)) and, for every l<h<r, MV{w{{Sj)j^h), A)) < My(((^j)jy/j), A("-''+i)). 
We take E such that M(X)M(wmax) < 

Finally, the complexity of Step 6 follows from ([5]). Note that, here, D < MV{A, A^""*^')) 
and L is the length of the sip computed in Step 4 according to equation (|4]) . □ 



4.3 Example 

Consider a sparse system of 2 polynomials in 5 variables supported on ^ = (^1,^2)) 

where Ai = {(0, 0, 0, 0, 0), (1, 1, 1, 0, 0), (2, 0, 0, 4, 2), (0, 0, 0, 8, 4)} and 

A2 = {(1)0, 1, 1, 2), (0, 1, 2, 5, 4), (1, 3, 0, 5, 4)}, with indeterminate coefhcients: 



f 



fi — Cu + C12X1X2X3 + CrsXfXfX^ + CuXfX^ 



5 

[/2 = C2lXiXsX4X^ + C22X2X^XlX^ + C23X1X2 

and the projection vr : — > A^, 7r(a;i, 2:2, ^3, X4, X5) = (xi,X2,X3). We are going to show 
a geometric resolution of 7r(y*(f)) following Algorithm /('-Projection. 

First, we apply Algorithm TransBasis and we obtain that {Xi, X2, X4} is a transcen- 
dence basis of K{V*{{)), and so, {Xi,X2} is a transcendence basis of K{Tr(V*{i))). 
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Now, the algorithm chooses a value b at random and specializes X4 = b. Set 6 = 1. 
The specialized system is 



fl : = 



fii = Cn + C12X1X2XS + CisXfX^ + C14XI 

/2I = (^21X1X3X1 + C22X2X^X^ + 



The next step is to apply Algorithm ParametricToricGeomRes with free variables 
Xi , X2 ■ We choose A = X^ as the primitive element to obtain the geometric resolution: 

•gjfsi-^i--^ +^^^ + c5 ^ + ^;cs ^ + 

2 



, -Cl2C2lCl3^!+C?lC22+Ci^2g23^?^l xA2 CliCjldlXf 
-\ 7T-7^ J- ^ 



•^^ ' Cl2^lX2 C12X2 Cl2^lX2 

w^iY) = Y 



Finally, Algorithm GeomResProj is applied to qx^jW^jW^ and a primitive element /u 
for K{7r{V*{f))). In this case, we take fi = X3 and obtain 

_ „ (v\ — v5 , __CrL_v4: I 2Ci2C23^i^|-Ci3C2iXf -^^3 2023(^22011X1+011014X1 -^^2 , 

• 9X3l^J-^ + O12X1X2 ^ + O12O22XI ^ + O12OI2XI ^ + 



_|_ Cl2C2gXfX2—Ci3C2lC2sXfY_^ OfgOllXiX, 
C12O22 *-^22*-^12 



4.4 Computation for generic rational coefficients. 

We are now going to show that, for a system with generic rational coefficients, the same 
steps as the ones in Algorithm i^-Projection can be performed using these rational 
coefficients to obtain a geometric resolution of the Zariski closure of the projection of the 
associated variety. 

Let {Xi, . . . , Xt,Xt+r+i, ■ ■ ■ ,Xn} be a transcendence basis of K{V*{f)) such that 
{Xi, . . . , Xf} is a maximal algebraically independent subset of {Xi, . . . , Xi}. 

As shown in [IH Appendix A], there exists a non-empty Zariski open set Ui C 
such that if c = (ci, . . . , c^) G ^1 n Q^, then: 

• the ideal Ic ■= (/i(c, X), . . . , /r-(c, X)) : (Xi . . . Xn)°° is radical equidimensional of 
dimension n — r, 

• {Xi, . . . , Xt, Xt+r+i, ■ ■ ■ , Xn} is algebraically independent modulo each of the asso- 
ciated primes of Ic, 

• for every 1 < k < i — t, {Xi, . . . , Xt, Xt^^} is algebraically dependent modulo Ic- 

Assume c Glii PlQ^ and consider the variety V*{i{c)) C A". Let W be an irreducible 
component of y*(f (c)). We have that dim(P^) = n — r and {Xi, . . . ,Xt, Xf+r+i, ■ ■ ■ , Xn} 
is a transcendence basis of Q{W). Then, the projection of W over the last n — t — r 
coordinates is a dominant map. Therefore, there is a Zariski open subset Ow C A""*"*^' 
such that, for every b G Ow n Q""*"'': 
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• Wb := wn{xt+r +1 — bt+r+1, ■ ■ ■ ,Xn — bn} IS ail equidimensional variety of dimension 
t, 

• {Xi, . . . ,Xt} is algebraically independent modulo I(Wb). 
Then, for every b e Ow n Q""*"'', the identity 7r{W) = 7r(Wb) holds. 

As the dimension of the set 9(f(c)) := V*{i{c)) - {x e (C*)'" : f(c,x) = 0}) is less 
than n — r, for every {ii, . . . ,it} C {1, . . . ,t + r}, there exists a non-zero polynomial 
Pii,...,it{Xii, ■ ■ ■ , Xi^., Xt+r+i, ■ ■ ■ 1 Xn) Vanishing identically on this set. Then, there is a 
non-empty Zariski open set Oi C A""*"*^' such that for every 6 € Oi, the dimension of 
0(f(c)) n {xt+r+i = bt+r+1, ■ ■ ■ ,Xn = bn} IS Icss than t. 

Then, for every 6 € Oi n fliy ^M/ H (Q*)"-*-^ we have that 

y*(f(c)) n {xt+r+1 = bt+r+1, ...,xn = bn} = {x G (C*)*+- : f (c, X, 6) = 0} X {6} 

and 7r(F*(f(c))) = [jw<W) = \Jw<Wb) = 7T{V*{i{c)) n {xt+r+l = h+r+l, ■ ■ ■ ,Xn = 
bn}); therefore, 

7r(y*(f (c))) = 7r({x G {C*Y+^ : f (c, f , 6) = 0} x {b}). 

For j = 1, . . . , r, let C (Z>o)*"'"'^' be the projection of Aj to the first t + r coordinates. 
Then fj = ^aes (^{aa)eA ^j,{a;a)-^'^)-^°^- Algorithm ParametricToricGeomRes works 
for generic sparse polynomial systems g supported on 5 = {Si, . . . ,Sr), that is, there 
is a polynomial pg in the coefficients of the system such that it computes a geometric 
resolution of V*{g{c)) for every vector of coefficients c with ps{c) ^ 0. Let U2 C 
be a non-empty Zariski open set such that, for every c := (ci, . . . , Cr), the polynomial 

i<j<r,aGSj) does not vanish identically. 

For c G n , there exists a non-empty Zariski open set O2 C A""*"^ such that 
for every b € O2 HQ"'"*"'', the algorithm ParametricToricGeomRes can be applied to the 
system f (c, x, 6). 

We conclude that, for coefficient vectors c G ^/i R ^2 H (Q*)^, a probabilistic algo- 
rithm Q-Projection which follows the same steps as Algorithm i^-Proj action can be 
applied in order to compute a geometric resolution of 7r(y*(f (c))). Taking into account 
the complexity estimates in Theorem [HI this proves Theorem [H 

Finally, we show an example where, following the steps of the algorithm Q-Projection, 
we obtain a geometric resolution of 7r{V*{f)) for a sparse system f with rational coefficients. 

Example. Let f be the following sparse system with the same support family A = 
{Ai,A2) as in the example of Section HTSl 

^ ^ f/i = 3 + 2X1X2X3 - XfXjXl + 5X|X| 

[f2 = 2X1X3X4X1 - 3X2X|X|X| + 7XiX|X|X4 

and let vr : — > be the projection it{xi,X2,x^,X4,X5) = {xi,X2,xs)- 

We use the same choices 6 = 1, A = X5 and /i = X3 as in Section HTSl and look at steps 

4 and 6 of the algorithm. 

In Step 4, the algorithm ParametricToricGeomRes computes the geometric resolution 

of V*{ii) with free variables Xi,X2 associated to the linear form A: 
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. qxs{y) = y 




~ 2x1X2 



2x1X2 



3 



• W5{Y) = Y. 



In Step 6, if Algorithm GeomResProj is applied to the geometric resolution obtained 
in Step 4 and the linear form ^, the geometric resolution i^Xz^v^i) is obtained, where 



This is, in fact, the geometric resolution of 7r(y*(f)) with free variables Xi, X2 associ- 
ated to the linear form /_i = X3, as it can be checked, for instance, by applying a Groebner 
basis elimination based procedure. 
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